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rom  four  of  the  early  nuclear  test  shots.  This  study  updates  the  parameters 
by  using  more  recent  and  thorough  experimental  data  and  by  using  information 
from  the  Department  of  Defense  Fallout  Prediction  System  (DELFIC  Computer 
Model) . 

Ps 

The  cloud  center's  height  is  updated  from: 

H (KFT)  = 44.0  + 6.1  InY  - .205  (InY  + 2.42)  I InY  + 2.421 
c 1 1 

to: 

Hc(KFT)  = 50.7  + 20.4  loglt)Y  + 3. 50(log1C)Y) 2 + 2.40(log1()Y)3  + 
0.60(log10Y)4 

The  Gaussian  spatial  activity  distribution  of  the  original  WSEG  model  is 
retained.  However,  both  the  vertical  and  horizontal  standard  deviations  are 
modified.  The  vertical  standard  deviation  is  changed  from  Ojj  = 0.180  Hc  to 
On  = 0.125  Hc.  The  horizontal  standard  deviation,  aQ,  is  multiplied  by  a yield 
dependent  adjustment  factor,  AK,  where 

AK  = 0.90  - 0.40  log1QY  + 0.  30(log1()Y) 2 + 0. 10(log1()Y) 3 

Inserting  the  modified  cloud  parameters  into  the  WSEG  model  results  in 
significantly  improved  fallout  predictions  below  500  kilotons  yield  when  DELFIC 
is  used  as  a standard  of  comparison.  Above  500  kilotons,  predictions  from  the 
updated  model  do  not  differ  significantly  from  the  original  WSEG  model. 
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Preface 


This  independent  study  began  as  an  effort  to  codify  the  basic 
nuclear  effects  on  the  TI-59  calculator.  Fallout  was  the  first  subject 
of  in-depth  study  with  the  WSEG/NAS  modified  fallout  model  selected  for 
codifying.  Critical  analysis  of  the  model  revealed  major  deficiencies 
in  fallout  prediction  capabilities.  The  deficiencies  appeared  to  be 
related  to  inaccurate  cloud  parameter  data  used  in  empirically  fitting 
the  original  model.  With  improved  cloud  information  available,  such  as 
from  the  DELFIC  prediction  model,  it  was  decided  to  refit  the  model's 
calculational  algorithms  and  compare  the  new  results  to  the  original 
model  as  well  as  to  the  DELFIC  model.  The  comparative  model  analysis 
confirmed  that  the  low  yield  capability  of  the  WSEG  model  is  significant- 
ly improved  by  the  corrected  geometry.  All  work  with  the  TI-59  project 
was  terminated. 

This  author  gratefully  appreciates  Professor  C.  J.  Bridgman,  thesis 
advisor,  for  his  catalytic  role  in  prompting  the  in-depth  study  of  fall- 
out which  resulted  in  these  finding.^.  Thanks  are  also  extended  to  Dr. 
Kathy  S.  Gant  and  the  Solar  and  Special  Studies  Section,  Oak  Ridge 
National  Laboratory,  for  providing  a copy  of  the  FORTRAN  version  of  the 
WSEG/NAS  model. 
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Abstract 


The  analysis  of  residual  radiation  from  the  detonation  of  nuclear 
weapons  has  resulted  in  a number  of  nuclear  fallout  prediction  models. 

The  most  widely  used  of  these  is  the  WSEG  (Weapons  System  Evaluation 
Group)  fallout  model  which  is  a series  of  empirically  based  analytical 
algorithms.  This  study  modifies  thpee  of  the  empirical  parameters  which 
control  the  transport  and  deposition  of  the  fallout  for  the  WSEG  model. 
These  three  are  the  height  of  the  center  of  the  nuclear  cloud,  the 
cloud's  vertical  thickness,  and  the  cloud's  horizontal  radius.  In  the 
original  model,  these  parameters  were  empirical  fits  to  data  from  four 
of  the  early  nuclear  test  shots.  This  study  updates  the  parameters  by 
using  more  recent  and  thorough  experimental  data  and  by  using  information 
from  the  Department  of  Defense  Fallout  Prediction  System  (DELFIC  Com- 
puter Model) . 

The  cloud  center's  height  is  updated  from: 

Hc(KFT)  = 44.0  + 6.1  InY  - .205  (InY  + 2.42)  | InY  + 2 . 42 j (1) 


H (KFT)  = 50.7  + 20.4  log  'f  + 3. 50(3 og. 0Y) 2 + 2. 40 (log, _Y) 3 + 

C J.U  J-U  xu 

0.60(log10Y)4  (2) 

The  Gaussian  spatial  activity  distribution  of  the  original  WSEG 
model  is  retained.  However,  both  the  vertical  and  horizontal  standard 
deviations  are  modified.  The.  vertical  standard  deviation  is  changed 
from  Oy  = 0.180  Hc  to  = 0.125  H . The  horizontal  standard  deviation, 
aQ,  is  multiplied  by  a yield  dependent  adjustment  factor,  AK,  where 

AK  = 0.90  - 0.40  log1QY  + 0.  30(logirY) 2 + 0. 10(log1()Y)  3 (3) 


Inserting  the  modified  cloud  parameters  into  the  WSEG  model  results 
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in  significantly  improved  fallout  predictions  below  500  kilotons  yield 
when  DELFIC  is  used  as  a standard  of  comparison.  Above  500  kilotons, 
predictions  from  the  updated  model  do  not  differ  significantly  from 
the  original  WSEG  model. 


AN  IMPROVEMENT  TO  THE  WSEG  FALLOUT 


MODEL  LOW  YIELD  PREDICTION  CAPABILITY 

I.  Introduction 

An  important  aspect  of  military  strategy  and  civil  defense  planning 
is  the  prediction  of  fallout  from  nuclear  detonations.  Current  limita- 
tions on  atmospheric  nuclear  testing  preclude  developing  predictions 
based  on  experimental  results  under  typical  battlefield  conditions. 
Therefore,  a number  of  empirically  based  fallout  models  have  been  devel- 
oped. The  most  widely  used  of  these  is  the  WSEG  fallout  model  which  is 
updated  in  this  study. 

As  early  as  1965,  it  was  recognized  that  the  WSEG  model  tends  to  be 
overly  pessimistic  in  predicting  casualties  from  residual  nuclear  radia- 
tion. Two  factors,  inaccurate  cloud  geometry  and  the  use  of  the  Way- 
-1. 2 

Wigner  T ' approximation  for  the  decay  of  the  fission  products,  were 
recognized  as  contributing  to  the  overprediction  of  fallout  by  WSEG 
(Ref.  1:9,207). 

More  recently,  Norment  has  suggested  that  the  deficiencies  of  the 
model  are  primarily  due  to  errors  in  the  activity  deposition  rate 
function,  g(t).  In  addition,  he  points  out  that  the  cloud's  altitude 
and  geometry  are  not  consistent  with  updated  data.  The  WSEG  nuclear 
cloud  is  too  low  and  too  small  in  radius  at  low  yields.  This  results 
in  fallout  contour  patterns  which  are  too  narrow  with  excessively  high 
concentrations  of  activity  too  close  to  ground  zero.  That  is,  the 
isodose  contour  lines  are  longer  and  narrower  than  those  observed  in 
test  shots  and  those  predicted  by  sophisticated  fallout  models  such  as 
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DELFIC.  It  is  in  the  low-yield  range  that  these  weaknesses  in  the  WSEG 
model  are  most  apparent.  (Ref.  2 : 84-90 ; 3: 29) 

This  study  has  applied  corrected  cloud  geometry  to  the  WSEG  model. 
The  corrected  cloud  is  higher  and  thinner  than  the  original  ,SEG  cloud. 
Also,  it  is  significantly  larger  in  radius  at  low  yields.  Since  the 
deposition  rate  function,  g(t),  is  dependent  on  cloud  geometry,  it  also 
has  been  improved.  The  improvements  made  to  the  model  are  validated  by 
comparison  to  the  DELFIC  model. 

The  comparisons  show  that  the  improved  model  is  much  more  consis- 
tent with  DELFIC  results  than  the  original  model  is.  This  is  a signifi- 
cant improvement  since  Norment,  in  a recent  evaluation  of  fallout  models, 
concluded  that  DELFIC  has  the  best  prediction  capability  for  six  test 
shots  considered  (Ref.  3:6). 

For  the  reader  not  familiar  with  the  parlance  of  fallout,  a list  of 
abbreviations  and  terms  common  to  fallout  literature  is  provided  in 


Appendix  A. 


II.  The  DELFIC  Model 


The  Department  of  Defense  Fallout  Prediction  System  (DELFIC)  is 
widely  recognized  as  the  most  credible  fallout  prediction  model.  How- 
ever, a widespread  lack  of  understanding  of  the  model,  combined  with  its 
requirements  for  large  amounts  of  computer  central  memory  and  computa- 
tion time,  tends  to  inhibit  its  use.  Appendix  B provides  supplementary 
information  on  DELFIC  for  the  reader  not  familiar  with  the  model. 

Norment,  in  a comparative  model  analysis  using  data  from  six  test 
shots,  recently  concluded  that  DELFIC  conforms  to  experimental  data 
better  than  any  other  prediction  model.  It  was  therefore  used  as  a 
comparison  standard  in  this  study.  Direct  test  data  was  not  used  in 
this  effort.  Norment' s evaluation  found  that  DELFIC  is  most  accurate 
up  to  50  kilotons  yield.  (Ref  3:6) 

A simplified  set  of  environmental  conditions  was  established  for 
obtaining  all  data  extracted  from  the  DELFIC  and  WSEG  models.  DELFIC 
was  used  with  Nevada  test  site  soil  and  with  the  U.S.  Standard  atmos- 
phere for  mid-latitude,  spring/fall  conditions.  A ground  roughness 
factor  of  0.5  was  used  for  all  accumulated  dose  comparisons.  A steady 
15  mph  wind  was  used  at  all  altitudes.  No  wind  shear  was  used.  The 
contour  width  plots  presented  give  the  normalized  unit  reference  dose 
rates  (Djj+^)  without  terrain  shielding  or  instrumentation  factors.  All 
contour  length  plots  use  the  ground  roughness  factor  of  0.5  and  give 
dose  accumulated  to  either  four  or  twenty  four  hours.  All  comparisons 
use  a 100%  fission  device  detonated  at  the  surface.  No  percentage 
figure  of  accuracy  is  assigned  to  the  three  models. 


III.  The  WSEG/NAS  Fallout  Model  (Ref.  2; 4; 5) 

The  version  of  the  WSEG  model  used  is  the  WSEG-10/NAS  modified  fall- 
out model  as  provided  to  Oak  Ridge  National  Laboratory  by  Leo  A.  Schmidt 
of  the  Institute  for  Defense  Analysis.  To  simplify  and  clarify  the 
nomenclature,  the  model  is  referred  to  simply  as  WSEG  throughout  this 
report.  WSEG,  in  its  various  versions  and  forms,  has  been  in  use  by 
analysts  for  many  years.  Although  its  results  can  be  in  error  due  to 
deficiencies  in  the  model,  systems  analysts  persist  in  using  the  model 
because  of  its  computational  ease  and  simplicity. 

The  WSEG  model  is  structured  around  an  empirically  based  deposition 
function,  called  g(t),  which  represents  the  rate  of  deposition  of  fall- 
out activity  on  the  ground  as  a function  of  time.  This  function  can  be 
conceptualized  in  the  following  way.  Consider  a nuclear  surface  burst 
which  results  in  a cloud  of  fallout  particles,  assumed  to  be  micro- 
spheres, at  an  altitude  well  above  the  ground.  The  center  of  this  cloud 

is  located  at  an  altitude  H . This  altitude  is  determined  by  the 

c 

weapon's  total  yield.  Consider  a no  wind  condition  so  that  the  cloud 
remains  stationary  over  the  burst  point.  The  following  four  paragraphs 
present  the  circumstances  and  processes  which  are  assumed  in  describing 
the  transport  and  deposition  of  the  fallout  by  the  WSEG  model. 

The  fission  product  radiation  and  unburned  fuel  were  fully  vapor- 
ized at  the  instant  of  burst  and  began  resolidification  into  fallout 
particles  as  the  cloud  rose.  The  resulting  solid  particles  will  have 
some  distribution  of  sizes  from  a few  microns  to  several  hundred  microns 
in  radius.  The  actual  distribution  depends  on  the  type  of  soil  which 
was  vaporized  and  on  the  weapon's  total  yield.  The  total  radioactivity 
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contained  in  and  on  these  fallout  particles  is  determined  primarily  by 
the  weapon's  fission  yield. 

Some  of  the  particles  will  fall  while  the  cloud  is  still  rising. 
This  is  known  as  stem  fallout  and  is  not  treated  by  the  WSEG  model.  The 
model  transports  and  deposits  local  fallout  which  is  roughly  80%  of  the 
activity  in  the  cloud.  The  remaining  20%,  consisting  of  the  smallest 
particles,  becomes  world-wide  fallout.  For  this  local  fallout,  the 
source  normalization  factor,  NF,  used  by  WSEG  is  2 x 10^  roentgens/hr/ 
fission  megatron/st  mi  normalized  to  one  hour  after  detonation.  This 
value  includes  both  fission  and  induced  activity.  For  bursts  above 
the  surface  of  the  earth,  the  activity  is  adjusted  by  the  fraction  of 
the  fireball's  lower  hemisphere  which  intersects  the  surface. 

The  various  size  particles  will  fall  from  the  stationary  cloud  with 
different  terminal  velocities,  raining  on  the  ground  below  at  varying 
times.  It  is  clear  that  there  is  some  function,  g(t),  which  describes 
the  deposition  rate  of  radioactive  material  on  the  ground.  This  rate  of 

material  deposition  does  not  include  the  effects  of  radioactive  decay. 

-1  2 

Decay  effects  are  applied  by  using  the  Way-Wigner  T ’ decay  rule. 

The  concept  of  a deposition  function,  g(t)  can  be  extended  to  the 
real  situation  with  wind.  When  wind  is  applied,  the  larger  particles 
reach  the  ground  somewhat  downwind  from  ground  zero  while  the  smaller 
particles  are  carried  to  greater  distances  downwind  from  ground  zero 
before  they  hit  the  ground.  Nonetheless,  the  function  g(t)  still 
describes  the  rate  of  arrival  of  activity  at  surface  level  as  a function 
of  time. 

It  is  clear  that  this  function  g(t)  is  dependent  on  the  height  of 
the  cloud,  Hc,  and  on  the  vertical  distribution  of  activity  about  Hc. 
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The  size  distribution  of  particles  and  the  activity  distribution  with 
various  sizes  also  influence  g(t) , as  do  the  density  and  viscosity  of 
the  atmosphere  between  the  cloud  and  the  ground.  Rather  than  attempting 
to  calculate  g(t)  from  all  these  factors,  WSEG,  using  four  of  the  early 
nuclear  test  shots,  approximates  g(t)  empirically  by  the  following 


equation. 


g(t)  = ~ exp  - (£  ) 


where  is  empirically  fit  as: 


Tp  = 1.0573203  Q _ 0<5  exp  _ (Hc)2-j  (5) 


where 


Hc  =44.0+6.1  InY  - .205  (InY  + 2.42)  |lnY  + 2.42 | 


(5. a.) 


The  g(t)  function  is  easily  transformed  to  a distance  function  g(l)  by 
applying  the  effective  wind,  EFW. 

g(l)  = ~ exp  - (1/L0)  (( 

"ft 


where 


Lq  = EEW.Tc  and  /“  g(l)  dl  = J*  g(t)  dt 


(6. a.) 


The  activity  within  the  nuclear  cloud  is  assumed  to  be  a Gaussian 
distribution  in  space.  This  distribution  is  supported  by  early  test 
data.  The  limits  of  the  distribution  are  defined  as  being  four  standard 
deviations  in  extent  both  vertically  and  horizontally.  The  density 
function  p describes  the  distribution  where  DWD  is  distance  from  ground 
zero  along  the  effective  wind  vector,  referred  to  as  the  hotline,  CWD  is 
distance  crosswind  from  the  hotline,  and  H is  activity  altitude. 

p(DWD, CWD,H)  = 7 T~  - y ~ (7) 

a ( °H  ) 


aG,  On»  and  Hc  determine  the  initial  dimensions  of  the  cloud.  The 
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original  empirical  fits  found  these  parameters  to  be: 


aQ(st.mi)  = exp  Qo.  70  + - 3. 25/(4. 0 + (InY  + 5.4)2)[]  (8) 

Hc(KFT)  = 44.0  + 6.1  InY  - .205  QlnY  + 2.42  ^ |lnY  + 2. 42  | (9) 

aH  - .18  Hc  (10) 

After  extensive  calculus  and  empirical  adjustments,  shown  in  part 
in  the  WSEG  report,  downwind  (f^)  and  crosswind  (f  ) transport  functions 
result. 


fd  = W*NF*F*g(l)*$  (Lo  ‘ DWD  } 

(L*V°d  ) 


(11) 


fc  = 1 exp  - \ ( QJD  ) 

<2ir)\.  (Vc  ) 


(12) 


W = Total  yield  in  kilotons 

2 

NF  = Source  normalization  factor  of  2000  r/hr/kt/st  mi 
F = Fission  fraction  of  the  device 

AF  = Fraction  of  the  fireball's  lower  hemisphere  intersecting  the 
earth's  surface 

g(l)  = the  deposition  function 

'J’(x)  = the  cumulative  normal  function  of  (x) 

L = (L  2 + 2a  2)^ 
o d 

2 2 2 2 2 2 
or  = a z (L  + 8 a )/(L  z + 2 a ) 
d o o o o o 

ax  = l/[l  + (.001  Hc  • EFW)/aQ2I 

°c2  ■ °o2  + r (8IDUD  + 2 °dl  °o2)  + h <0dTc°usc)2  + r* 

((DWD  + 2 od)LoIcOHSc)2 

Sc  = Wind  shear  component 


This  model  uses  a wind  shear  over  the  vertical  extent  of  the  cloud. 
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The  shear  component  (Sc)  is  the  change  in  the  wind  component  normal  to 
the  effective  wind  within  the  cloud  divided  by  the  cloud  thickness.  It 
is  entered  in  mph/kft  of  cloud  thickness.  No  downwind  shear  is  used  in 
the  model. 

The  normalized  unit  reference  dose  rate,  in  roentgens  per  hour,  at 
the  point  (DWD,CWD),  is  the  product  of  the  downwind  and  crosswind 
transport  functions. 

Dh+1(DWD,CWD)  = fc*fc  (13) 

-1  2 

The  T * decay  approximation  is  applied  to  the  activity  in  its  transit 
downwind. 

An  effective  biological  dose  is  also  determined.  This  dose,  which 
is  assigned  a multitude  of  names  in  the  literature,  will  be  referred  to 
as  simply  the  biological  dose.  The  dose  calculation  assumes  that  ten 
percent  of  the  dose  received  is  not  repairable  and  that  the  remaining 
ninety  percent  is  repairable  with  the  damage  decaying  over  a thirty  day 
time  constant. 

The  biological  dose  is  approximated  from  the  following  expression. 

Bio  Dose  ~exp  - {.287  + .52  In  g + .4475  In  ^6)2}  (14) 

Ta  is  the  average  fallout  arrival  time.  The  expression  is  obtained  by 
integrating  the  unit  reference  dose  rate  from  the  time  of  arrival. 

D(t)  = 10%  f Dh+1  x'1-2dT  + 90%  /tDH+1  t"1-2  exp 
3 Ta 

{(30  DAYS) (t-l) } dx  (15) 

Four  Fortran  subroutines  for  the  WSEG  fallout  model  are  reproduced 
in  Appendix  C.  No  main  program  or  provisions  for  constructing  contours 
are  provided  since  each  user's  requirements  are  unique.  No  conversion 


to  the  metric  system  is  made  and  the  subroutines  are  presented  in  their 
received  form.  The  model's  outputs  are  stored  in  an  array  as  follows: 
ARRY(33)  is  biological  dose;  ARRY(34)  through  ARRY(39)  relate  to  special 


military  requirements;  ARRY(23)  is  fallout  arrival  time;  ARRY(32)  is  D . 
Additional  user  instructions  appear  as  comments  in  each  subroutine. 

In  this  study,  the  three  empirically  determined  parameters  are  re- 
computed on  the  basis  of  more  complete  experimental  data  and  on  DELFIC 
calculations  to  update  the  information  from  the  original  four  early  test 
shots.  These  three  are:  Hc,  the  height  of  the  cloud’s  center;  0H,  the 
vertical  standard  deviation  of  the  activity  distribution  in  the  cloud; 
and  0O,  the  initial  value  of  both  0^  and  0C  which  determine  the  hori- 
zontal extent  of  the  cloud. 

New  values  for  these  three  parameters  are  developed  in  Part  IV. 

The  effects  of  these  new  values  on  the  WSEG  model  are  presented  in  the 
comparative  model  analysis  of  Part  V. 
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V 


Previous  evaluations  of  fallout  models  have  recognized  that  the 
WSEG  nuclear  cloud's  height  and  dimensions  are  inconsistent  with  infor- 
mation which  is  more  complete  and  current  than  the  four  early  test  shots 
used  to  empirically  fit  WSEG  (Ref  1:93;  3:29).  In  addition,  Norment  has 
recently  proposed  that  the  deposition  rate  function,  g(t),  requires  up- 
dating (Ref  2:86).  In  the  original  WSEG-10  report  published  in  1959, 
Pugh  and  Galiano  stated  that  their  empirical  values  are  based  on  limited 
data  and  that  they  should  be  revised  when  improved  information  becomes 
available  (Ref  5:17).  However,  except  for  a minor  adjustment  to  the 
empirical  Tc  in  a WSEG-10  supplement,  this  author  could  locate  no 
correction  or  refinement  to  the  values  assigned  to  the  empirical  para- 
meters in  the  original  report.  Therefore,  this  study  updates  the  cloud 
parameters  by  using  more  current  information. 

Four  different  sources  of  information  for  the  cloud  parameters 
were  evaluated  for  use.  Recognizing  the  merits  and  limitations  of  each 
source  made  it  necessary  to  quantitatively  evaluate  each  parameter  from 
each  source  independently. 

The  first  set  of  cloud  parameters  chosen  for  evaluation  were  the 
values  of  the  original  WSEG  model.  Although  the  cloud  was  known  to  be 
too  low  and  dimensionally  in  error,  it  was  selected  as  a basis  of  com- 
parison for  all  parametric  changes. 

The  second  source  of  cloud  information  was  data  from  a series  of 
reports  following  completion  of  all  atmospheric  testing.  These  reports 
were  prepared  by  the  old  Defense  Atomic  Support  Agency  (DASA).  Many  of 
the  cloud  parameter  curves  in  these  reports  are  simple  empirically 
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based  extrapolations  using  power  functions  of  yield.  Much  of  the  data 
is  based  on  visual  observations.  All  data  used  from  these  reports  was 
extracted  from  Ref  2 by  Norment.  For  simplicity,  all  information  from 
these  reports  will  be  called  DASA  data. 

The  third  source  of  cloud  information  used  was  Effects  of  Nuclear 
Weapons,  ENW  (Ref  6:431).  The  ENW  cloud  parameters  result  from  an 
empirically  based  atmospheric  model. 

The  DELFIC  model  was  the  final  source  evaluated  for  use.  The 
recently  revised  Atmospheric  Sciences  Associates'  version  of  DELFIC  was 
used  for  all  DELFIC  information  appearing  in  this  study.  All  dose  and 
dose  rate  predictions  by  the  original  and  updated  WSEG  models  are  com- 
pared to  DELFIC  predictions. 

In  selecting  the  best  parametric  fits,  the  most  emphasis  was  placed 
on  the  range  of  accumulated  dose  from  approximately  fifty  rads  to  about 
two  thousand  rads.  This  range  is  important  to  the  user  studying  the 
multiple  burst  scenario  as  well  as  to  the  civil  defense  planner  involved 
in  shelter  studies.  The  stem  fallout  predictions  of  DELFIC  were  ignored. 
Prediction  curves  in  the  region  of  lethal  prompt  effects  were  also 
ignored.  All  study  was  limited  to  the  range  of  yield  from  one  kiloton 
to  thirty  megatons.  The  low-yield  region  of  highest  DELFIC  credibility 
was  studied  most  closely. 

The  Cloud  Thickness  Correction 

The  original  WSEG  model  uses  a cloud  thickness  of  40y  where  CTy  = 
0.180  Hc.  Norment  has  recently  recommended  that  a = 0.125  Hc  would 
be  more  consistent  with  experimental  data  from  DASA  (Ref  2:84).  Cloud 
thickness,  using  both  0.125  Hc  and  0.180  Hc,  was  plotted  versus  yield 
for  the  ENW  and  WSEG  models'  cloud  heights.  DELFIC  thickness  from  the 
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cloud  rise  module  was  also  plotted  on  these  same  graphs  which  appear  as 
Figures  1 and  2.  Observe  the  exceptional  agreement,  using  .125  Hc,  for 
all  three  curves  at  the  lower  yields.  The  original  WSEG  value  of  0.180 
Hc  is  not  consistent  with  either  DELFIC  or  ENW  models.  Although  the 
0.125  Hc  value  is  recognized  as  an  improvement  to  the  original  value  of 
.180  Hc,  both  values  will  be  evaluated  in  correcting  the  other  parameters. 
The  Cloud  Height  Correction 

In  empirically  determining  the  height  of  the  cloud  center,  Pugh 
and  Galiano  used  the  base  of  the  visual  cloud  as  the  center  of  the  active 
cloud  (Ref  5:24).  This  is  now  known  to  be  too  low.  The  top  of  the  WSEG 
cloud  is  1.36  Hc,  using  Ojj  = 0.180  Hc,  and  1.25  Hc  using  the  improved 
Ojj  = 0.125  Hc.  These  are  plotted  in  Figures  3 and  4 respectively.  Also 


l. 

plotted  are  the  DELFIC,  DASA,  and  ENW  cloud  tops.  Several  test  shots' 

• j] 

visible  tops  are  also  included  (Ref  1:85).  Since  these  shots  in  the  IVY 
series  were  under  atypical  conditions,  with  coral  and  water  drawn  into 
the  cloud,  they  are  presented  merely  for  informative  purposes  (Ref  1:85). 

Note  the  tendency  for  the  visually  based  DASA  cloud  top  to  be 
excessively  high  at  most  yields.  Up  to  about  five  megatons,  the  DELFIC 
and  ENW  models'  cloud  tops  agree  exceptionally  well.  Above  that, 

DELFIC' s cloud  top  exceeds  all  other  predicted  and  observed  tops.  This 
excessive  cloud  height  is  consistent  with  Norment's  finding  that  DELFIC 
underpredicts  activity  at  the  higher  yields  (Ref  3:5).  This  divergence 
of  DELFIC' s cloud  top,  as  well  as  its  thickness,  at  the  higher  yields  is 
due  to  the  DELFIC  cloud  rise  module's  treatment  of  the  vertical  extent 
of  the  cloud  as  a point  within  the  atmosphere.  The  known  tendency  of 
the  WSEG  cloud  to  be  low  is  most  apparent  only  at  the  higher  yields  when 
the  top,  rather  than  the  center,  of  the  clouds  are  compared. 
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CLOUD  THICKNESS  VS  YIELD 
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CLOUD  THICKNESS  VS  YIELD 
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Figure  2.  Cloud  Thickness  vs  Yield  for  CJH  = .125  H 
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Cloud  Top  vs  Yield  for  aH  = . 18  H 


CLOUD  TO?  VS  YIELD 


Based  on  the  preceeding  discussions  and  on  Figures  3 and  4,  it 

initially  appeared  that  the  ENW  model  cloud  top  provides  the  best  cloud 

top  information.  Before  this  was  adopted  for  use,  however,  the  DASA, 

ENW,  and  original  WSEG  cloud  tops  were  tested  in  the  model.  They  were 

tested  using  both  ct,  of  .180  H and  o of  .125  H . 

H c H c 

The  for  the  DASA  and  ENW  cloud  tops  were  determined  by  setting 

the  cloud  top  value  for  each  equal  to  1.36  H for  a„  = . 18  H and  equal 

c H c ^ 

to  1. 25  H for  aTI  = . 125  H . All  four  H values  were  computed  and 
c H c c r 

plotted  versus  yield.  Then,  the  four  were  least  squares  fit  to  calcu- 
lational  algorithms.  The  resultant  algorithms  are: 

DASA  DATA,  a„  = .18  H 
H c 

Hc(KFT)  =53.0  + 19.1  log1()Y  + 3. 40(log10Y) 2 + 1.  40(log1()Y)  3 + 

.30(log10Y)4  (16) 

DASA  DATA,  au  = .125  II 
H c 

Hc(KFT)  = 57.  7 + 20.8  log^Y  + 3.  70(log1()Y)2  + 1. 40(log1()Y) 3 + 

.30(log1()Y)4  (17) 

ENW  MODEL,  a„  = . 180  H 
H c 

Hc(KFT)  = 46.6  + 18.8  log1()Y  + 3.20(log1()Y)2  + 2. 20(log1{)Y) 3 + 

.60(log1()Y)4  (18) 

ENW  MODEL,  0H  = .125  Hc 

Hc(KFT)  = 50.7  + 20.4  log1()Y  +*50  (log^Y) 2 + 2. 40(log1()Y)  3 + 

.60(log10Y)4  (19) 

Each  of  the  four  H values  were  then  inserted  into  the  WSEG  model  with 
c 

H* 


the  appropriate  cr. 


Accumulated  doses  were  calculated  using  the  DASA  data  and  ENW  model 
Hc  values  with  of  .125  Hc  and  .180  Hc.  The  algorithms  were  tested 

. ] 

from  one  kiloton  to  thirty  megatons  of  yield  and  compared  to  DELFIC  pre- 

( dictions.  Typical  results  are  seen  in  Figures  5 and  6.  The  original 

WSEG  model  is  also  plotted  in  those  figures.  It  was  observed  that  as 
the  yield  increased  through  500  kilotons,  the  effects  of  the  changes  to 
vertical  geometry  resulted  in  minimal  changes  to  predictions. 

Throughout  the  yield  range  considered,  the  DASA  Hc  algorithm  tended 
to  underpredict  close-in  activity  and  overpredict  activity  further  down- 
wind. This  is  consistent  with  its  excessive  cloud  height  in  Figures  3 
and  4.  This  finding  supports  Russell’s  early  work  which  concluded  that 
visual  cloud  observations  do  not  necessarily  reflect  active  cloud  para- 
meters (Ref  1;93).  Therefore,  the  DASA  information  was  rejected  for  the 
new  Hc  determination. 

The  original  WSEG  and  DELFIC  model  cloud  tops  were  rejected  for  use 

. 

because  of  the  arguments  presented  previously. 

Throughout  the  yield  range,  it  was  seen  that  the  ENW  model's  Hc 
with  Oy  = .125  Hc  consistently  predicted  fallout  contour  lengths  which 
were  shorter,  thus  conforming  better  to  DELFIC  predictions,  than  did  the 
ENW  model  with  = .18  Hc.  Thus,  the  initial  conclusion,  that  using 
the  ENW  model  cloud  top  with  a CJjj  = .125  Hc  is  the  best  available  cor- 
rection to  the  vertical  cloud  parameters,  was  confirmed  as  correct. 

The  updated  algorithm  for  the  height  of  cloud  center  is: 

Hc(KFT)  = 50.7  + 20.4(log10Y)  + 3.50(log10Y) 2 +2.  40(log1C)Y)  3 + 

0.60(log10Y)4  (20) 

The  corrected  = .125  Hc.  This  corrected  vertical  geometry  was 
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incorporated  into  the  WSEG  model. 

The  Cloud  Radius  Correction 

Since  Russell  has  shown  that  the  visual  cloud  radius  observed  in 
test  shots  does  not  accurately  describe  the  active  radius  of  the  cloud, 
BELFIC  radius  was  used  exclusively  in  updating  the  WSEG  model's  radius 
(Ref  1:93).  Figure  7 shows  stabilized  cloud  radius  versus  yield  for 
DELFIC  and  for  WSEG.  The  DELFIC  radius  plotted  is  the  radius  immediately 
following  stabilization  of  the  top  in  the  cloud  rise  module.  Considera- 
tion was  given  to  using  the  laterally  stabilized  DELFIC  radius,  but  that 
radius  resulted  in  insignificant  changes  to  predictions  at  and  below  100 
kilotons.  Above  100  kilotoris,  the  fallout  contours  were  grossly  widened 
and  reduced  in  length.  Those  wider  and  shorter  contours  are  more  con- 
sistent with  DELFIC  results  at  the  higher  yields,  but  DELFIC  has  been 
validated  only  up  to  fifty  kilotons  and  is  known  to  under-predict  at  the 
higher  yields.  In  addition,  it  has  been  found  that  WSEG  radius  is 
excessively  small  only  at  lower  yields  and  is  adequate  at  high  yields. 
Figure  7 shows  that  the  vertically  stabilized  DELFIC  cloud  radius  used 
is  in  good  agreement  with  that  result. 

From  one  kiloton  to  thirty  megatons,  a multiplying  factor  for  WSEG 
radius  was  calculated.  Multiplying  WSEG  radius  by  the  factor  results  in 
the  DELFIC  radius  value.  This  multiplying  factor,  called  AK,  was  then 
plotted  versus  yield  and  also  appears  in  Figure  7.  AK  was  fit  by  least 
squares  into  a yield  dependent  algorithm. 

AK  = .90  - .40  log1QY  + .30(log1()Y)2  + .10(log1C)Y)3  (21) 

Multiplying  the  aQ  in  the  WSEG  model  by  AK  corrects  WSEG  radius  to  agree 
with  DELFIC  radius.  Use  of  the  AK  multiplying  factor  results  in  a maxi- 
mum deviation  of  corrected  WSEG  radius  from  DELFIC  radius  of  16%  in  the 


l ' 

. i 

1 


21 


STABILIZED  RADIUS  VS  YIELD 
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two  megaton  region.  Higher  and  lower  yield  regions  remain  well  below 
that  deviation. 

The  data  studied  in  this  section  and  the  comparative  model  analysis 
of  the  following  section  show  that  the  WSEG  model  is  more  sensitive  to 
changes  in  horizontal  geometry,  or  radius,  rather  than  to  changes  in 
vertical  geometry. 

The  three  empirical  changes  made  to  the  WSEG  model  in  this  section 
were  incorporated  into  the  Fortran  version  of  the  model  in  Appendix  C. 
The  revised  Fortran  subroutine  FALLY  is  presented  in  Appendix  D. 


V.  Comparative  Model  Analysis 


The  original  WSEG  and  corrected  WSEG  models  are  both  compared  to 
DELFIC.  The  twenty  four  hour,  and  some  four  hour,  accumulated  doses 
along  the  downwind  hotline  are  plotted  versus  distance.  Typical  low 
and  high  yield  results  are  seen  in  Figures  8 and  9.  Additional  data 
appears  in  Appendix  E.  Note  how  the  WSEG  model's  excessively  long  iso- 
dose contours  are  shortened  significantly  at  low  yields.  For  example, 
at  5 kilotons,  the  100  rad  24  hour  contour  line  is  reduced  from  35  miles 
downwind  for  the  WSEG  model  to  only  20  miles  for  the  updated  model. 

This  is  close  to  the  DELFIC  length  of  16  miles.  When  each  model's 
errors  are  considered,  the  DELFIC  and  updated  WSEG  curves  are  in  very 
good  agreement  outside  of  the  stem  fallout  region.  This  good  correlation 
between  the  updated  model  and  DELFIC  is  observed  through  the  range  of 
maximum  DELFIC  validity,  one  through  fifty  kilotons.  Above  that,  the 
updated  and  original  WSEG  models  both  exceed  DELFIC  predictions.  Note 
that  in  the  100  to  200  kiloton  yield  range,  the  updated  model  lies 
between  the  DELFIC  and  WSEG  predictions.  Qualitatively,  this  compares 
favorably  to  the  results  of  test  shot  Koon.  WSEG  overpredicts  Koon  and 
DELFIC  underpredicts  it  (Ref  3:5).  In  general,  DELFIC  tends  to  under- 
predict at  the  higher  yields.  Figures  1 and  5 support  this  conclusion 
when  the  effects  of  excessive  cloud  height  and  thickness  are  considered. 

At  and  above  500  kilotons,  the  predictions  of  the  corrected  and 
original  WSEG  models  are  virtually  indistinguishable  within  the  bounds 
of  model  accuracy  up  to  about  10  megatons.  Note  the  effect  of  the 
differences  in  radius  between  the  original  and  improved  models  in  the 
one  to  ten  megaton  range.  These  radius  effects  cause  the  predictions 
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of  the  updated  model  to  exceed  those  of  the  original  model  only  in  this 
area.  The  sensitivity  of  the  model  to  radius  changes  can  be  seen  by 
comparing  Figure  7 to  the  Appendix  E data  and  the  crosswind  comparison 
charts  of  Appendix  F,  Above  10  megatons,  the  updated  model  again  tends 
to  shorten  and  widen  the  fallout  contours  consistent  with  Figure  7. 

The  effects  of  the  updated  cloud  geometry  on  contour  width  are  seen 
in  Figures  10  and  11.  Additional  data  is  given  in  Appendix  F.  Note,  in 
Appendix  F,  that  the  updated  model  consistently  predicts  wider  contours 
than  the  original  model  does  at  all  yields  where  WSEG  has  been  critiqued 
as  having  contour  lines  which  are  too  narrow. 

At  low  yields,  where  DELFIC  has  been  shown  to  be  most  valid  and 
where  WSEG  has  been  shown  to  be  most  unreliable,  the  shortened  and  wid- 
ened contours  of  the  corrected  model  more  closely  resemble  DELFIC  pre- 
dictions, which  have  been  shown  to  best  model  actual  experimental  data, 
than  the  contours  of  the  original  WSEG  model  do.  (Ref  3) 

Above  50  kilotons,  quantitative  assessment  of  the  improved  model  is 
not  feasible  due  to  the  lack  of  test  data  under  such  conditions  and  the 
lack  of  an  experimentally  validated  comparison  model  at  the  higher 
yields.  The  divergence  of  the  updated  model's  radius,  and  therefore 
contour  widths,  from  DELFIC  and  original  WSEG  in  the  two  to  ten  megaton 
range  is  therefore  considered  as  inconsequential  when  compared  to  the 
confirmed  improvements  made  at  the  low  yields  where  the  improvement  can 
be  quantified  and  validated. 
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Figure  10.  7 Kt  Isodose  Width  Comparison 
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VI.  Conclusions  and  Recommendations 


By  correcting  the  nuclear  cloud's  dimensions  and  height,  the  WSEG 
fallout  model's  prediction  capability  at  yields  below  500  kilotons  has 
been  significantly  improved.  In  this  age  of  radically  improved  inertial 
guidance  systems,  MIRV  for  ICBM  and  SLBM,  the  extremely  accurate  cruise 
missile,  and  reduced  emphasis  on  the  manned  penetrating  bomber,  it  can 
be  expected  that  the  strategic  planner  will  place  future  emphasis  on 


highly  accurate,  lower  yield  weapons  for  optimum  utilization  of  his 
available  throw  weight  and  weapon  material.  Thus,  the  improvements  made 
to  the  WSEG  model  are.  highly  relevant  to  the  yield  range  of  most  prac- 
tical interest.  Therefore,  the  effects  of  the  geometrical  changes  and 
their  calculational  algorithms  to  WSEG  in  the  multi-megaton  range,  which 
cannot  be  quantitatively  evaluated  by  either  typical  experimental  data 
or  by  a verified  prediction  model,  will  be  considered  as  being  of  merely 
pedagogical  interest. 

The  corrected  model's  improved  prediction  capability  at  the  low 
dose  levels,  below  500  rads  accumulated,  is  highly  applicable  to  the 


scenario  of  many  lew  yield,  highly  accurate  weapons  detonated  in  a 
target  area  rather  than  one  or  a few  high  yield  weapons.  Any  over- 


prediction, as  is  done  by  the  original  WSEG  model,  in  the  total  (from 
all  weapons)  dose  accumulated  range  of  300  to  600  rads  can  grossly 
exaggerate  casualty  predictions.  Very  small,  relative,  overpredictions 
in  the  LD/50  region  around  450  rads  accumulated  can  cause  disproportion- 
ate overprediction  of  casualties.  Russell  has  proposed  that  WSEG,  in  an 
earlier  version,  could  conceivably  overpredict  fatalities  by  one  to  two 


This  study  recommends  that  users  of  the  WSEG  fallout  prediction 


model  make  three  changes  to  it.  The  height  of  the  cloud  center,  Hc, 
should  be  corrected.  The  current  model's  cloud  center  height  is  given 
by: 

H (KFT)  = 44  + 6.1  InY  - .205(lnY  + 2.42)  IlnY  + 2.421  (22) 

c 

This  should  be  changed  to: 

Hc(KFT)  = 50.7  + 20.4  log1()Y  + 3. 50(log10Y) 2 + 2.  40(log1()Y)  3 + 

.60(log10Y)4  (23) 

This  correction  is  seen  as  ARRY(3)  in  Appendix  D,  the  corrected  sub- 
routine FALLY.  The  cloud  thickness  correction  is  made  by  changing  a 

ri 

from  a value  of  .180  H to  .125  H . This  OtT  is  ARRY(4)  in  Appendix  D. 

c c H 

The  radius  correction  is  made  by  multiplying  the  0Q,  or  ARRY(l),  by  a 
yield  dependent  adjustment  factor,  AK  in  Appendix  D. 

AK  = .90  - .40  log1QY  + .30(log1()Y)2  + .10(log1()Y)3  (24) 

By  making  these  three  changes  to  the  original  WSEG  model,  the  WSEG 
user  can  significantly  improve  prediction  capability  below  500  kilotons. 
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Appendix  A 


List  of  Terms  and  Abbreviations 


This  appendix  provides  the  reader  who  is  not  familiar  with  the 
nomenclature  used  in  this  study  with  a list  of  common  terms  and 
abbreviations  used  herein. 


AFIT 

DELFIC 


ENW 

I CBM 

IVY 

LD/50 


MIRV 

NAS 

RAD 


Roentgen 

SLBM 

WSEG 


Air  Force  Institute  of  Technology 

Department  of  Defense  Land  Fallout 
Prediction  System 

Unit  Reference  Dose  Rate.  Gives 
activity  at  detonation  plus  one  hour 
counting  all  fallout  that  has  arrived 
at  and  is  in  transit  to  the  point  of 
interest. 

EFFECTS  OF  NUCLEAR  WEAPONS,  widely  used 
reference  by  S.  Glasstone  and  P.  Dolan, 
see  reference  6. 

Intercontinental  Ballistic  Missile 

Test  shot  series  on  Pacific  atolls 

Exposure  level  in  rads  defined  as  the 
level  where  essentially  100%  fatalities 
result  at  twice  LD/50  and  very  few 
occur  at  one  half  the  LD/50  dose. 

Multiple  Independently  Targeted  Re- 
entry Vehicle 

National  Academy  of  Science 

Radiation  absorbed  of  100  ergs  of 
ionizing  radiation  per  gram  of  absorb- 
ing material 

Dated  term  approximately  equivalent  to 
one  rad 

Submarine  Launched  Ballistic  Missile 
Weapons  Systems  Evaluation  Group 


Appendix  B 

The  DELFIC  Model  (Ref  4:9-10) 


This  appendix  is  presented  to  provide  the  reader  unfamiliar  with 
the  DELFIC  fallout  prediction  model  with  some  general  information  about 
the  model.  Since  the  model  is  far  too  detailed  and  complex  to  explain 
indepth,  the  model  summary  by  H.  Norment,  in  Reference  4,  is  given  in 
its  entirety. 

"DELFIC  is  a research  code  which,  for  practical  purposes,  is  useful 
only  to  those  who  have  the  time  and  inclination  to  become  deeply  involved 
in  local  fallout  prediction.  It  is  structured  in  a physically  straight- 
forward manner  such  as  to  include,  via  use  of  the  best  practicable 
models,  all  of  the  phenomena  that  are  important  to  the  formation  and 
distribution  of  local  fallout  from  surface  and  low  airbursts.  It  is 
designed  for  highly  flexible  usage. 

The  code  uses  a dynamic  cloud  rise  model  that  produces  results 
which  are  demonstrably  superior  to  those  produced  by  conventional  means. 
This  model  has  the  unique  advantage  of  being  able  to  account  for  the 
effects  of  atmosphere  structure  on  the  cloud  rise  and  stabilization. 

Thus,  it  can  made  credible  predictions  for  locations  that  are  geographi- 
cally remote  from  the  test  site  areas,  where  all  of  the  observed  data 
have  been  obtained.  (Examples  of  important  remote  areas  are  northern 
Europe  and  northern  Asia.)  Activity  calculations  are  rigorously  done 
such  that  use  of  a questionable  activity  normalization  factor  (i.e., 

the  K-factor)  is  not  required.  Also  not  required  is  use  of  a gross- 

-1. 2 

activity  decay  equation,  such  as  the  conventional  t ' law. 

Transport  can  be  via  time  and/or  space  variant  wind  fields,  or  it 


14 


can  follow  the  practice  of  all  other  codes  by  using  a single  wind  field 
that  varies  only  in  the  vertical.  Turbulent  dispersion  is  described  in 
Appendix  A. 4. 

Maps  of  a large  number  of  different  properties  of  the  fallout  field 
can  be  prepared,  including  some  types  that  are  not  available  from  any 
other  code. 

In  ref.  3 we  describe  two  versions  of  the  code:  one  used  by  the 
Ballistics  Research  Laboratories  (BRL) , and  an  updated,  improved  version 
used  by  Atmospheric  Science  Associates  (ASA).  In  the  ref.  3 study  the 
BRL  version  was  used,  but  here  we  use  the  ASA  version. 

In  the  ref.  3 study  we  found  that  DELFIC  provides  adequate  pre- 
dictions for  surface  and  near  surface  bursts,  and  we  rated  it  to  be  the 
best  of  the  codes  studied." 
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Original  WSEG  Model  Fortran  Subroutines 


I NE  FBLLY  74/74  0PT  = 1 FTN  4.D+446  09/13 


SUBROUTINE  FALL Y <YI EL C,FISS, HOB, ARRY) 

C TO  COMPUTE  THE  yiflT  DEPENDENT  °A  ? AM  ET  EPS  IN  THE  WSEG/0/-NAS 

C FALLCUT  MODFL.  YIELD  IS  Y T FLO  IY  MEGATONS , FISS  IS  FISSION 

C FRACTION,  H 03  IS  HEIGH-  OF  3 UP ST  IN  FEET,  ARRY  IS  AN  AWRY  OF  FORTY 

C ELEMFNTS  USED  TO  ORE3FPVE  t?FSUL-5  OF  DIFFERENT  SUBROUTINE  CALLS. 

C SUQRGUTINES  SHOULD  3E  fjlLEO  IN  ORDER  OF  NEW  YIEIO,  WIND  VELOCITY 

C OR  WIND  SHEAR,  DOWNWIND  DISTANCE,  AYD  CFO SS HIND  DISTANCE. 

C THE  VALUES  IN  AR^Y  < 1 ) TO  ARRY (6)  ARE  FILLED  HERE. 

DIMENSION  ARPY(l) 

XLNY  = ALOG (YIELD) 

TEM=XLNY+5.4 

TEMP= 0.70+0. 333  3 S3  3«  X LN  Y -3 . 25/ < 4 . 0 + T EM*  T EM ) 

ARRY  (1)  = EXP(TEMP) 

ARRY ( 2 ) = ARRY(l)  ♦ARPY(l) 

TEMP  = XLNY  +2.4? 

• C ORIGINAL  W'EG/NAS/S/°  VEPSION 

ARRY  (3)  = 44.+  =.i*XLNY  -0 . 20 5“ T EM?* APS ( T EMP) 

ARRY ( 4 )=. 13+ARRY (3) 

HCTVO  = ARRY ( 3)  /?f>. 

HCSIX  = ARRY ( 3 ) / : 0 • 

ARRY <5 >=1.057 320? • ( 12.* HCSIX-2 ,c * HD 5 1 X«HCSI X > * ( 1 . 0- 0. 5»EXP (-HCT WO* 
1HCTWO) 

IF (HOB . GT . 0 . ) GO  TO  ft 
ARRY(f)  i 2000000. *cISS*YIELO 
RETUF  N 

6 CONTINUE 

XMHB=160.» (YIELD* 1000. )*+0.4 
IF (HCB.LE.XMH9)  30  TO  10 
ARRY ( 6) =0 • 

RETUF  N 
10  CONTINUE 

TEHP=HOB/XMH  P 

AF=0.5« (1 .-TEMP)  + ( 1 .-TEMP) M 2.+TEMP)  + 0.001* TFMP 

ARRY (6)=200 000  0 .’  FITS* AF* YIELD 

RETURN 

FNO 


o o o o o 


* 


> SUPEOUTIN'  FALL'WIEFW.SC.ARRY) 

TO  COMPUTE  THE  W I »1 0 SPrEO  DP  HIND  SHEAR  EPPECTS  IS  THE  HSEG  10- 
NAS  FALLOUT  MODEL. 

EFW  IS  WIN)  V'LJDI TY  If:  STATUTE  MILES  PER  «OUR,  SC  IS  WIND  SHEAR 
IN  S'ATUIE  MILES  3PR  HOUR  PEP  KTLD'OOT,  AR'Y  IS  A STORAGE  ARRAY. 
THE  VALUES  IN  ARRY CM  to  ApPY  < 1 3)  ARE  SUPPLIED  HERE. 

PI MEN  SION  ARRY(l) 

XLO  = Er  We  AkR Y ( ? ) 

XL  OS  = XLO»  XLO 

SIGUS  = ARRY ( ? ) * (X.OS*3.’A  RRY ( 2) )/(<LDS+2.*ARRY (?) ) 

ARPY(IA)  = SORT  (SI  G'JS) 

XLS  = XLCS  ♦ 2. ‘SIGNS 
XL  = SORT  (XLS) 

AP.PY  < 15)  = 6./X. 

ARRY  ( 7 ) =APRY  (1  !>)  • A RRY  ( ?) 

TMPA  r ARP.Y(C)4A  = RY(A)'SC 
TEMP  = XLOMMPA/XL3 
ARRYO)  = TEIP'-TEWP 

TEMP  = ARRY (ID  • ARRY (C)*ARRY(A> * SD/XL 
ARRY(S)  = A RRY ( 2)  TE HP* T EMP 

XLOOS  = XLOS  ♦ 0.5*3IG"S 
ARRY  (10)  = (XLO'  ♦ SK-IJS)  /XLOPS 
IF(A,RY(10)  .LT. 1,002)  GO  TO  8 
TM  = 1 ./ARRY (10) 

C GAMMA  FUNCTION  tpc^ox  HASTINGS  ’.156 

GAMMA  = 1.  HMM-0. 57869867  ♦ TM* ( 0 . 9? 7 6 17 8 1 *• TH* ( - 0 . 6 23S£>27*-TM* ( 
10.67399080  + T M*  ( - 0 . 3?.?  2793  * TM»  0 . 0 7 67 320 £ ) ) ) > > 
ARRY(11)=1./(XL*SAMNA) 

GO  TO  9 

8 ARRY f 1 1 ) * 1 . /XL 

9 CONTINUE 

ARRY ( 17)  = 0.001'APRY(T)*EFW/ARRT(1) 

' ALONF  = 1./  (1.  «■  ARRY  ( 17)  ) 

«’  ARRY  (18)  = XLO/  ( X L*  AL  ONE*  ARP.Y  ( TUI) 

ARRY  (12)  = XLOS  * A RP  Y ( 5 ) * A RRY  ( 5 ) / ( X LS*  XL  O PS ) 

ARRY (13)  = 0.25  ♦ 2. ' SIGUS/XLOPS 
IE(EFW.LT .0.000001)  GO  TO  5 
ARRY ( 16)  = 2./EcW 
GO  TO  6 

G ARRY ( 10)s9999999i3. 

6 CONTINUE 
RETUFN 
END 


******* 
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6 

7 


c 


?1 


SUBROUTINE  FALLOW) 040,  HOC ALjTWSN, ARRY) 

TO  COMPUTE  OOWNWIND  DISTANCE  EFFECTS  IN  THE  WSEG  10  - NAS  FALLOUT 
MODEL . 

OKO  IS  TWF  OOWNWINO  DISTANCE  IN  STATUTE  MILES,  MDCAL  = O-ONLY 
WSEG  BIO  DOSE,  1-NMCSSC  °I0  DOSE  AT  BMP 3 , 2-ALSO  MAY  DOSE, 

TWPN  IS  TIME  OF  WEAPON  DETONATION  IN  HOURS,  ARRY  IS  A STORAGE 
ARRAY . 

THE  VALUES  IN  AR?Y(?0)  TO  ARRY (31)  ARE  COMPUTED  HERE. 

DIMENSION  ARRY(l) ,RRRS(5> 

DATA  PHPS  ( 1 ) , °HR  S ( ?)  , nHRS ( 3 ) , HRSC4),  BHF.S  (5)  / 

1 7.,??.,  hi.,  211 .,  800./ 

TP  = 0WD*2.*A-!RY  ( 1 !•  > 

DWP=ABS(TP> 

TMP  = ARRY ( 15)  *DWP 
IE<TMP.LE.3 .)  GO  TO  8 
0WP=3 • /ARRY ( 15) 

8 CONTINUE 

SIGCS  = ARRY  (9)  «•  ARRV(7)  *OWP  ♦ ARP.  Y <6 ) * TP*  TP 
SIGC  = SORT (SIGCS) 

TA  = ARRY ( 16) * OWD 
IF(TA.GT.4.)  GO  TD  11 
APPROY  HASTINGS  P.185  FPR  CUM  NO" 

TM  = ABS(TA/1  .41121356'’) 

TMP  = 1.  ♦ TM* (0.278393+TM* ( 0 . 230339 +TM* (0. 000972*TM*0. 078108) ) ) 

TMP  = TMP  ‘TMP 

CUP  = l.-l./fTM"' T M° ) 

IF (TA  .LT.  0)  GO  TO  5 
CUV  = 0.5*(1.  * CUP) 

GO  TO  7 
CONTINUE 

CUV  = 0.5* ( 1 . - CUP) 

CONTINUE 

ALTWC’  = l./(l.+APRY(17)*(l.-CUV)  I 
GO  TO  12 

11  ALTHO= 1. 

1 2 CONTINUE 

ARRY <?0>=l./<  2.50553* SIGC) 

TMP  = ALTWO  » SI GC 
ARRY ( 21 ) = 0.5/(TMP«TMr) 

TA  = OWD  * ARRY ( 18 ) 

IFUA.LT.5.)  GO  TD  14 
CUV=1 . 

GO  TO  15 
14  CONTINUE 

APPROY  HASTINGS  P.187  FOR  CUM  NOR. 

TM  = APS ( T A /I  .411.21356?) 

„ TMP  = 1+TM*  ( 0.  07P52307D4*TM»  (0.  01223  201234-TM*  (0. 0 092705272  + TM* 
1(0.0001520143*TM' ( 0 . 0 0 02765672*71 » 0 . 0 0 00* 3 Ofc 38 ) ) ) ) ) 

TMP=TMP*TMP 
TMP=T  MF*TMP 
TMP=TMP*TMO 


tmp=tmf*tmp 

CUP=  l.-l./TMP 

IF  (TA  .LT.  0)  GO  TO  21 

CUV  = 0.5* ( 1 . ♦ CUP) 

GO  TO  22 
CONTINUE 


iO  XO  DOC 
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CUV  = 0.5M1.  - CUP) 

22  CONTINUE 
15  CONTINUE 

IF(APRY(10) .LT.l.OO?)  GO  TD  17 
TMP  = <ABS(0WD)*»REY(i;;)/8.)**ARRM10) 

GO  TO  18 

17  TMP=ABS(DW9)*ARRY(15)/6. 

18  CONTINUE 

IF (TUP. LT .30.)  GO  TO  19 
ARRY (22)  = 0. 

ARRY (23)  = 999995. 

ARRY  ( 24 ) = 1. 

RETURN 

19  CONTINUE 

GT  =ARRY(11)* EX°(-TMP) 

ARRY  (22)  = ARRY(i)  •9T*f'UV 
TMP  = ARRY ( 1 3)  * TP  * T P*ARRY ( 1?) 

ARRY (23)  = SQRT(T“P) 

TMP  = ALOG(AR<vY(23>/31.6) 

ARRY (24)  = EXP( -(0.  23 7 ♦ 0 . 52*TMP»0. 0947 5*TMP* TMD> ) 

IE (MDCAL  .NE.  0)  GO  TO  31 
PETURN 

31  CONTINUE 

TMP  = ARRY(23)**(-0.2) 

DO  32  J = 1,5 
JST  = 24  ♦ J 
BT  = BHRS(J)  - TWPN 
BTT  = BT  - ARRY(?3> 

IF  (BTT  .GE.  0.)  GO  TO  33 
ARP.Y  ( JST)  = 0. 

GO  TO  32 

33  CONTINUE 

77  = 0.5  + 4.5*EXP(-(  0.00051  ♦ 0. 0 00 25*T HP) * ( BTT > ) 

ARRY ( JST)  = (TMP  - (3T«* (-0.2)) •*?7 
3?  CONTINUE 

IE (MPCAL  .NE.  1)  50  TO  34 
RETURN 

34  CONTINUE 
TMP  = 0. 

00  35  K = 1,5 
KLK  = 24  ♦ K 

IF ( TMP  .GT.  ARRY ( KLK) ) GO  TO  35 
TMP  = ARRY(KLK) 

- KVL  = KLK 

35  CONTINUE 
IF (KVL  .NE.  23)  50  TO  36 
ARRY (30)  = ARR Y ( 25 ) 

ARRY ( 31 ) = BHRS(l) 

GO  TO  39 

36  CONTINUE 
IF (KVL  .NE.  29)  50  To  37 
ARRY ( 30)  = APRY(?9> 

ARRY (31)  = BHRS ( 5 ) 

GO  TO  39 

37  CONTINUE 
YM  = ARRY (KVL  - 1) 

YO  = ARRY (KVL) 

YP  = ARRY (KVL  ♦ 1) 

OELLT  = 0.25* (AL0G(PHRS(5) ) - AL05 ( BHRF ( 1) ) ) 

TO  = A L 0 G ( B H RS (KVL  - 24  ) ) 

DELP  = YP  - YO 
HELM  = YO  - YM 
OELSO  = DELP  - OFLM 
DEL YC  = 0 . 5*  ( DEL  B"  ♦ DELM) 

OT  = - DELYO*OELLT/OELSQ 
XL  T = TO  ♦ 

ARRY ( 3 1 ) = 

ARRY ( 30)  = 

33  CONTINUE 
RETURN 
FNO 

THIS  PAGE  IS  BEST  vtu«***- - — 
39  FROM  COP  I FURNISHED  10  DDC 


OT 

EXP ( XL  T ) 

YO  - 0 •5*0ELY0*QFLY0OF.S'I 


«„M.tTYPRACTlCABia 


OOOnOOOODOOOO 


SUBROUTINE  FALLCH(CWD,»OCAL , APRY) 

TO  CCMFUT  E CROSS'UNT  0 TST  AMC  E EF-EOTS 
FALLOUT  MODEL  AND  ORODUCE  FINAL  AHSHE 
CWD  IS  CPOSSWINO  01  STANCE  TN  ST'TJTE 
910  TOSE,  1-NMCSSO  time  ooses,  2-ALSO 
ARRAY. 

FOP  OUTPUT  the  H + 1 OOSE  RATE  IS  in 
PICLCGICAL  OOSE  IS  IN  ARRY (33),  HE  T 
AFTEF  WEAPON  DURST  TIMF  IS  IN  ARRY  < 2 3 
THE  NMCSSC  810  OOSE  AFTER  7,22, E3, 211 
STRAT  OF  THE  TIME  AXIS  IS  IN  ARRY ( 3H ) 
PIOLCGICAL  POSc  IS  IN  ARRY(39>  AMO  TH 
TIME  ORIGIN  IS  IN  AR»Y<31>. 

THE  VALUES  IN  ARRY ( 3?)  TO  ARRY(33)  AR 
DIMENSION  ARRY(l) 

TMP  = CWD  * A RRY (21) * C WO 
IF  (TRP  .GT.  30.)  GO  TO  6 
FC  = ARRY (20)  *EX°(-Trf°> 

ARPY  (32)  = FG *APRY (22) 

CONTINUE 

ARRY ( 33) = ARRY (32) *ARRY(24> 

IF ( MOCAL  .NE.  0)  GO  TO  3 
RETUFN 
CONTINUE 
00  (*  J = 

ARRY ( J ♦ 

CONTINUE 
IF (MOCAL  .NE. 

RETURN 
CONTINUE 

ARRY (39)  = ARRY(30)*ARRY (32) 

RETURN 
CONTINUE 
ARRY ( 32) 

GO  TO  7 
FNO 


FOR  THE  WS EG  10  -NAS 
RS. 

MILES,  MOCAL  OF  O-ONLY  WSEG 
MAX  DOSE,  ARRY  IS  A STORAGE 

ARRY (32) , THE  WSEG 
IMF  OF  FALLOUT  ARRIVAL 

) .- 

, AND  800  HOURS  FROM  THE 
TO  ARRY (38) . THE  HAX 
E TIME  OF  MAX  OOSE  AFTER 

E COMPUTED  HERE. 


1,5 

33)  = A5RY(32)»ARRY (J  t- 


1)  GO  TO  5 


?'+> 


= 0. 


eP 


& 


y 
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Appendix  D 

Improved  Model  Fortran  Subroutine  FALLY 


SUBROUTINE  FALL Y( YIELD, FI SS. HOB, ARRY) 

TO  COMPUTE  THE  YIELO.DEPENCi.HT  .PARAMETERS  IN..THE  K3EG/0/-N„S  

FALLOUT  MOOEL*  Y I ELC  13  YIELD  IN  MEGATONS,  l'ISS  IS  FISSION 
FRACTION,  HOB  IS  HEIGHT  OF  BURST. IN  FEET,  ARRY  IS  AN  ARRY.  OF  FORTY 
ELEMENTS  USED  TO  PRESERVE  RESULTS  C-  DIFFER  EE  if  SUBROUTINE  CALLS. 

SUBROUTINES  .SHOULD  BE  CALLED.  IN._OKlER._CF  MlW  YlttD,  Wl  M..VELO.CXXX 

OR  HIND  SHEAR;  OOWNWINO  DISTANCE,  AN'S  CR0S3WIN0  DISTANCE. 

THE  VALUES  IN  ARRY  ( 1)  TO... ARRY ( fcJ.  AE£_EILL£D  HERE. 

DIMENSION  ARRY ( 1 ) 

XLNY  = ALOG  (YIELQ1 1 

TEM  = XLNY  *b . A 

AK=  . 9-.  A*  ALOG L0.(_YIRLQ1.+  .3*  LA  LOG  10  ( YIE  LO)  ) *»2.  * ,lMALQG.lfl.iYJI£LDJJ 

1**3, 

TEMP  = 0.?  0+G. 333333  i*  XL  NY.- 3.  2 5 /.  I h.t.  0 *-LEH*.T  E tf  J_ 

ARFY  (11  =•  EXP(TEMP) 

ARRY  ( 1)  = AK*  ARRYY  1 ) _ 

ARFY (2)  = ARRY(l)  *ARRY ( 1 ) 

TEMP  =.  X LNY_T.2j.-42- 

C IMRROVEO  ENW  FIT 

ARF  > 13)  = 50. 7+20.V*ALOGL0  ! Y.IE LD ) T3  . 5*_I_AL3G1 0 (Y.I ELQ).).f  *2  . + 2.  4.* 

l(ALOG10CYIELU))*->3.t.6*(AL0G10(YIELD))**A. 

C. SIGH  A H IS  .12  5. -OF.-H  _C 

5 ARK Y ( 4)  = • 1 2 5* ARRY ( 3) 

HCTWO  =_.AF_R5ti.31  E2S. 

rtCSIX  = AR}>  Yl3)  /w(J» 

A RR Y (E ) si. 0 5732  3 3* ( 12 . *HCSIX -2 . E» HCS1 X*HCS1X) * (1 . 0- 0 . 5 * EXP ( -HCTHO* 
1HCTW0)1 

IF  ( H03.GT.  0 «.)  GO  _TO.  6 

ARRY (6)  = 2000000. *FISS»YIELO 

RETURN 

6 CONTINUE 

XMHB=  130  .♦  ( YIELD*  1 000  . ) **0«.4 - 

IF (HOB.LE. XMHB)  GO  TO  10 

ARRY  (6)  - 0. .• : 

RETURN 

— 10.  CONTINUE 

TEMP=H03/XMHB 

AF=0.5*(  l.-TEMP)  * ( 1. -TEMP)  *(■*.  + TEMP)  + 0.0 01* TEMP 

ARRY (61 =2000003. *FISS*AF» YIELD 

..  RETURN 
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Appendix  E 


Isodose  Contour  Length  Comparisons 

The  graphs  on  the  following  pages  plot  accumulated  dose  in  rads 
versus  distance  downwind  for  a number  of  yields.  The  original  WSEG 
model,  the  updated  WSEG  model,  and  DELFIC  predictions  are  presented  for 
comparison.  Isodose  contour  lengths  are  easily  compared  for  the  three 
models  by  selecting  an  accumulated  dose  level  on  the  vertical  axis  and 
then  reading  horizontally  to  the  right  until  each  of  the  three  curves 
is  intercepted.  The  mileage  value  on  the  horizontal  axis  directly 
below  the  intercept  point  is  the  maximum  length  of  the  contour  for  that 
model.  All  data  is  for  100%  fission  yield  devices  detonated  at  the 
surface  of  the  earth.  The  ground  roughness  factor  is  .5. 
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Appendix  F 

Isodose  Contour  Width  Comparison 

This  appendix  presents  contour  width  comparisons  between  the 
original  and  improved  WSEG  models  and  DELFIC  for  several  yields  at 
selected  downwind  distances.  No  ground  roughness  factor  is  applied. 
All  data  is  for  100%  fission  yield  devices  detonated  at  the  surface  of 
the  earth.  The  output  plotted  is  the  normalized  unit  reference  dose 
rate  (D^) . 
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